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Abstract. We propose a nonlocal theory of single-particle excitations. It is based on an off- 
diagonal effective medium and the projection operator method for treating the retarded Green 
function. The theory determines the nonlocal effective medium matrix elements by requiring 
that they are consistent with those of the self-energy of the Green function. This arrows 
for a description of long-range intersite correlations with high resolution in momentum space. 
Numerical study for the half-filled Hubbard model on the simple cubic lattice demonstrates that 
the theory is applicable to the strong correlation regime as well as the intermediate regime of 
Coulomb interaction strength. Furthermore the results show that nonlocal excitations cause sub- 
bands in the strong Coulomb interaction regime due to strong antiferromagnetic correlations, 
decrease the quasi-particle peak on the Fermi level with increasing Coulomb interaction, and 
shift the critical Coulomb interaction Uc2 for the divergence of effective mass towards higher 
energies at least by a factor of two as compared with that in the single-site approximation. 



The description of single-particle excitations in correlated electron systems has been a 
challenging problem because new phenomena are observed there and because simple perturbation 
approaches are not applicable [Ij. At present, excitations in the systems of infinite dimensions 
have been clarified by use of the dynamical mean field theory (DMFT) [2]. That theory makes use 
of a momentum-independent self-energy, and determines the energy self-consistently by solving 
an impurity problem with use of advanced many-body theories. The DMFT can be traced back 
to the many-body coherent potential approximation (CPA) for magnetic alloys [3]. It is also 
known to be equivalent to the dynamical CPA [4] for describing finite-temperature magnetism 
[5]. More recently, we proposed the projection operator method CPA (PM-CPA) |6], which is 
equivalent to the DMFT too. It is based on the projection technique applied to the retarded 
Green function and on an energy-dependent Liouville operator. 

The theories mentioned above are based on a single-site approximation (SSA) which neglects 
the off-diagonal components of the self-energy. Therefore the momentum dependence of 
excitations in real systems is described incompletely. In this paper, we propose a nonlocal theory 
of excitations which extends the PM-CPA, and demonstrate that it describes self-consistently 
momentum-dependent excitations with high resolution in the strong Coulomb interaction regime 
as well as in the intermediate regime. 

We adopt here the Hubbard Hamiltonian H with nearest-neighbor transfer integral t and 
intra-atomic Coulomb interaction U, and consider the retarded Green function 



G.,.(^) = (al|-^«y . (1) 
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Here the inner product between two operators A and B is defined as (^|-B) = ([A''", with 
use of a thermal average ( ) and the anti-commutator [ , ]+. Furthermore z = uj + i6 with 6 
being an infinitesimal positive number, and L is a Liouville operator defined by LA = [H,A]^ 
for an operator A. 

According to the Dyson equation, the Green function is expressed as 

Gijjz) = [{z -Ho- A{z))-%^ . (2) 

Here {Hojiju is the Hartree-Fock Hamiltonian matrix and {K{z))ija- is the self-energy matrix 
defined by Aijfj{z) = U'^Gijaiz). The reduced memory function Gija{z) is given by 

G,,„{z) = [aI I • (3) 

The operator A\^ is defined by A\^ = aj^Sni^rr with 5nj_o- = nj_o- — {ni^a-)- (ojo-) is the 
creation (annihilation) operator for an electron with spin a on site i, and {riia) is the average 
electron number on site i for spin a. Moreover, L is a Liouville operator acting on the space 
orthogonal to the space {|at)}; L = QLQ, Q = 1 - P, and P is a projection operator defined 

by -P = Eia |aL) (alal- 

In our nonlocal theory, we introduce an energy-dependent Liouville operator L{z) for 
an effective Hamiltonian Hq{z) of an off-diagonal medium «-e., Hq{z) = Hq + 

X^jjCT ^ija{z) aj^aja- Here Hq is the Hartree-Fock Hamiltonian. The retarded Green function 
Fija{z) for the Liouvillean L{z) is expressed as 

F,,,{z) = [[z-Ho-nz))-%, , (4) 

where {T,{z))ij(j = Tjijfj{z). It should be noted that the Green function Fij„(z) becomes identical 
with Gija{z) when 

^ija = ^ijo-{z) • (5) 

Thus, Fiju{z) describes properly many-body excitations when the above relation is satisfied. 

In order to obtain an explicit expression of the self-energy Aiju{z), we separate the Liouvillean 
L into L{z) and the remaining interaction part Li{z), i.e., L = L{z) + Li(z), and expand the 
memory function ([3]) by using the incremental method: 

= G^{z)+T.^f^^i^+\T. E AG^iz) + .... (6) 

Here Acg? (z) = G^!^ (z) -Cfji! (z), and Acfff^ (z) = -^f^ (z) - Acfff (z) - Aof^ (z) - 

Gij^{z). These terms are all calculated from cluster memory functions defined by Gij^{z) = 

{A\^ |(z A]^). l}-''\z) = QL^''\z)Q and L(^)(z) (c = ...) is the Liouvillean for 

a cluster c embedded in the medium with off-diagonal matrix elements {S^mo-C-z)}- Note that a 
"cluster" c = (ij), for example, does not mean that sites are nearest neighbors. Instead 
they may be far apart. 

When we truncate the expansion at a certain stage in the series of Eq. ([6]), the self-energy 
Ajjo-(-z) = U'^Gija{z) depends on the medium Tiij^{z). We then determine the medium self- 
consistently from condition ([5]). Note that the present theory reduces to the PM-CPA (i.e., the 
DMFT) in infinite dimensions. 



Figure 1. Momentum-dependent spectra along high-symmetry hnes at U = 14. The Fermi 
level is indicated by a bold line. 



In the above self-consistent theory, one needs explicit expressions for the cluster memory 
functions, for which we apply the renormalized perturbation theory (RPT) [7]. 



Gg) (z) = iG'^'iz) ■ (1 - Lr ■ G'^'{z))-'\,, . (7) 



Here {Lj = 11(1 - 2{ni^^))/xia, and Xia = {ni-a){l - {ni-a))- In the simplest 

approximation (RPT-0), {Gq )ija{z) is given by 

1*^0 )ij(T{Z) — ^ija / _ / // ■ (oj 

Here = Xia/ {ni^a)c{l - {n^-a)c)Sij + 1 - Si,, x{e,e',e") = (1 - /(e))(l - f{e'))f{e") + 
/(e)/(e')(l — /(e")), /(e) is the Fermi distribution function, (njo-)c is the electron number for a 

cavity state defined by {m^), = ^ depfl{e)f{e), and ^^(e) = -^-^lTn[{F,{z)-^ + t!-'\z))-%^. 

The coherent cluster Green function {Fc{z))ija = Fijfj{z) is given by Eq. ([J]), and (Xl {z))ija = 
T,ij„{z) for sites belonging to the cluster c. 

Finally the momentum dependent excitation spectra are calculated from the Green function 

GUz) = \ f , ■ (9) 

Here Cka is the Hartree-Fock one-electron energy eigen value, and the momentum-dependent 
self-energy is calculated via Fourier transform of the off-diagonal self-energy as Akf^{z) = 
J2j ^jOa{z) exp(ifc ■ Rj). 

We have performed numerical calculations for the single-particle excitation spectra on a 
simple cubic lattice at half- filling. In the nonlocal self-energy calculations, we have taken all 
two-site pairs up to 10-th nearest neighbors, and neglected the clusters beyond the 2 sites. Figure 
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Figure 2. Densities of states in the present 
theory (soHd curve) and the SSA (dotted 



curve 



1 shows an example of the momentum-dependent spectra along high-symmetry lines for U = 14 
in unit of \t\ = 1. We find quasiparticle states near the Fermi level, which form a narrow band 
with averaged quasiparticle weight Z = 0.21. On the other hand, the lower (upper) Hubbard 
band is enhanced around the T (R) point. It should be noted that the energy splitting between 
the main peaks of the lower and upper Hubbard bands is about 20, being larger than U = 14 
as seen in Fig. 2. This seems to be explained by strong antiferromagnetic (AF) correlations. In 
fact, in the strong AF correlation limit the energy to remove (add) an electron is expected to 
be eo — zJ {eo + U + zJ) at half filling. Here eo is the atomic level and J is the super exchange 
interaction constant (J = 4\t'\^ /U). Then the splitting is given hy U + 2zJ instead of U. The 
former yields 17.4 for U = 14, which is comparable to the calculated splitting 20. We also find 
sub-peaks around |a;| = 8 in Figs. 1 and 2. These peaks are interpreted as the lower and upper 
Hubbard peaks in the absence of long-range AF correlations because the splitting is close to U . 
Furthermore, we find weak excitations at ~ 3, which might correspond to a "shadow band" 
due to strong AF correlations as found in QMC calculations [8j. We also point out that in the 
nonlocal theory, the density of states at the Fermi level p(0) does not agree with the one in the 
SSA as seen in Fig. 2. 

We have calculated the nonlocal excitation spectra with increasing Coulomb interaction. Self- 
consistent solutions are found up to C/ ~ 21. In this range, we find that the quasiparticle weights 
Z are larger than those of the SSA but never vanish. We suggests that the critical Coulomb 
interaction is Uc2{Z = 0)^30. 

In summary, we have presented a nonlocal theory for excitation spectra by introducing a 
fully off-diagonal effective medium Tiij(j{z) to the projection operator formalism. The theory 
describes self-consistently the long-range nonlocal excitations in momentum space with high 
resolution. We verified from the numerical calculations that the present scheme works from the 
weak to strong Coulomb interaction regime, and found the shadow bands as well as the Hubbard 
sub-bands due to strong AF correlations. 
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